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Heavy Flavour Working Group are shown, using four different theoretical 
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1 Introduction 



Measurements of inclusive charmless semileptonic decays of B mesons, B Xuii^, 
are directly related to the CKM matrix element \Vub\- The theoretical description [1] 
of the hadronic current involved in these decays, relying on the Operator Product Ex- 
pansion (OPE) technique, allows the determination of \Vub\ from the total decay rate 
with a small uncertainty. However, in order to suppress background from semilep- 
tonic decays with charm, B — Xciiy, measurements of partial branching fraction are 
performed in restricted kinematic regions. Unfortunately, OPE breaks down in some 
of these regions, and the theoretical uncertainty increases significantly. Contributions 
due to weak annihilation also play a role in part of the kinematic regions, and need 
to be carefully assessed. In short, theory and backgound subtraction give conflicting 
requirements, and a trade-off must be found. Due to the improved knowledge of 
B — )• Xcii' transitions and to the abundant data samples collected at the B factories, 
recent results based on phase space regions which are increasingly larger allow for an 
improved precision in the \Vub\ determinations. 

Experimental measurements of charmless semileptonic decays are reviewed in Sec- 
tion [2], with emphasis on new preliminary results from Babar. Preliminary \Vub\ 
averages from the Heavy Flavour Averaging Group by using the available theory 
calculations are presented in Section |3i Conclusions are given in Section HI 

2 Measurements of Partial Branching Fractions 

Measurements in the endpoint region of the lepton momentum spectrum are concep- 
tually simple, being based on the identification of an high momentum lepton only. 
However, the kinematic region selected by the high lepton momentum requirement 
to suppress charmed background suffers from sizeable theoretical uncertainties. All 
experimental efforts [2j have been devoted to reducing the lepton momentum cut as 
low as allowed by background knowledge. Signal-to-background ratios (S/B) of the 
order of 1/10 have been achieved, with signal efficiencies at the 30% level or less. 
An improved analysis [3], based on the measurement of missing energy to estimate 
the maximum kinematically allowed hadronic mass squared, s™"^, resulted in S/B of 
about 1/2. A summary of the available endpoint measurements is given in Table [1] 
The recoil technique aims at fully reconstructing one of the two B mesons (Breco) 
from the T(4S') decay in a fully hadronic decay, which allows to determine completely 
the decay kinematics of the other B (Brecon)- It is therefore possible to access relevant 
kinematic variables, such as the invariant mass of the hadronic system, mx, the light- 
cone momentum component P+ = Ex — \px\, and the squared invariant mass of 
the lepton pair, g^. Semileptonic events are identified by an high-momentum lepton 
{p} > 1 GeV) and a missing mass consistent with zero. Non semileptonic backgrounds 
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Experiment 




E, (GeV) 


Ai3(10-^) 


El 


Babar 


81.4 


2.0-2.6 


5.72 ±0.41 ±0.65 


E2 


Belle 


27.0 


1.9-2.6 


8.5 ± 0.4 ± 1.5 


E3 


CLEO 


9.13 


2.2-2.6 


2.30 ±0.15 ±0.35 


E4 


Babar (^r^^ <3.5 GeV^) 


81.4 


2.0-2.6 


4.41 ±0.42 ±0.42 



Table 1: Summary of endpoint analyses, labeled E1-E4 in the following. 



Kinematic Region 


Signal Yield 


AB{B XJv) (10-3) 


Babar Rl 


< 1.55 GeV/c 


1033 ± 73 


1.08 ±0.08 ±0.06 


R2 


Mx < 1.70GeV/c 


1089 ± 82 


1.15 ±0.10 ±0.08 


R3 


P+ < 0.66 GeV 


902 ± 80 


0.98 ±0.09 ±0.08 


R4 


R2 and > 8 GeV/c^ 
p} > IGeV/c, (Mx,g2) fit 


665 ± 53 


0.68 ±0.06 ±0.04 


R5 


1441 ± 102 


1.80 ±0.13 ±0.15 


R6 


p} > 1.3GeV/c 


562 ± 55 


0.76 ±0.08 ±0.07 


Belle R5 


p} > IGeV/c, (Mx,g2) fit 


1032 ± 91 


1.96 ±0.17 ±0.16 



Table 2: Summary of signal yields and partial branching fractions in six kinematic 
regions (labeled R1-R6 in the following). Unless otherwise noted, the lepton momen- 
tum in the center-of-mass frame is required to be p| > 1 GeV/ c. The uncertainty on 
the yields is statistical only. 



are subtracted by studying the distribution of the beam-energy substituted mass 
for the -Breco candidates. About 1000 hadronic modes are reconstructed, with 
efficiencies at the 0.3% (0.5%) level for neutral (charged) B decays. 

Background from B — )■ XcCv is reduced mainly by vetoing charged kaons and Ks-, 
whose production is highly suppressed in signal, and charged and neutral soft pions 
kinematically compatible with B — )■ D*lv decays. 

Both B Factory experiments published results using the recoil technique [3] . Babar 
recently released [5J a preliminary update based on the full dataset (426 fb~^), which 
is detailed in the following. The integrated luminosity analysed by Belle is 626 fb~^. 

Partial rates for charmless semileptonic decays have been measured in several 
phase space regions, defined in Table as well as for charged and neutral B decays 
separately. The latter is achieved by explicitly requiring the (absolute) Breco charge 
to be one or zero, respectively, after subtacting with Monte Carlo a small fraction 
of events where the Breco charge was not correctly reconstructed. In addition to the 
mx, P+ and {mx,q^) distributions, the lepton momentum p*^ was also studied. A 
background-enriched control sample, obtained by reversing the vetoes on kaons and 
soft pions mentioned above, was used to determine the relative contribution due to 
semileptonic decays into P-wave D mesons directly on data; Monte-Carlo simulation 
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Figure 1: Upper row: Mx (a), P+ (b), with Mx < 1.7 GeV/c^ (c) and p| (d) 
spectra (data points), measured in Babar data. The result of the fit to the sum of 
three MC contributions is shown in the histograms: B — > X„£z/ decays generated 
inside (no shading) and outside (dark shading) the selected kinematic region, and 
B — > Xc^v and other background (light shading). Lower row: corresponding spectra 
for B — )■ Xudv (not corrected for efficiency) after background subtraction. 



was then reweighted accordingly. Although the impact on signal yields was almost 
negligible, the fit chisquares improved significantly. 

The event yields and partial branching fractions are given in Table [2l The distribu- 
tions of the kinematic variables under study, before and after background subtraction, 
are shown in Figure [H Statistical uncertainties range from 7% to 9%. 

A summary of the systematic uncertainties is given in Table [3l which also shows 
the corresponding uncertainties from the Belle recoil analysis. The statistical and 
experimental systematic uncertainties are of the same order. Detector-related uncer- 
tainties are dominated by undetected or mismeasured particles {e.g. Ki and addi- 
tional neutrinos) from background. Progress on the knowledge of exclusive B — )■ X(.lv 
decays refiects in a relatively small uncertainty due to background composition. In the 
most inclusive phase space region (R5), the dominant uncertainty is due to the signal 
model, in particular to the knowledge of heavy quark parameters and the branching 
fraction of exclusive B — )■ X^^v decays, which are used in the simulation to determine 
signal efficiency. Total uncertainties range between 9% and 13%. 

Measurements of partial rates for charged and neutral B mesons allow to determine 
the relative contribution of weak annihilation to the total rate, 7vi/A/r. The resulting 
90% confidence level regions are reported in Table H] they are in agreement with 
previous determinations |6]. 
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Source 


Rl 


R2 


Babar 
R3 R4 


R5 


R6 


Belle 
R5 


Statistical error 


7.1 


8.9 


8.9 


8.0 


7.1 


8.9 


8.8 


MC statistics 


1.3 


1.3 


1.3 


1.6 


1.1 


1.2 




Detector-related: 


2.8 


3.7 


5.5 


4.1 


3.2 


2.7 


3.3 


Fit-related: 


2.7 


4.9 


3.2 


3.2 


2.1 


2.5 


3.6 


Signal model: 


2.7 


3.0 


3.5 


1.9 


6.6 


7.9 


6.3 


Background model: 


2.0 


2.6 


3.4 


2.8 


2.8 


2.2 


1.7 


Total systematics: 


5.3 
-5.0 


6.4 
-6.2 


8.0 
-8.1 


6.2 

-6.2 


8.5 
-7.7 


9.4 

-8.7 


±8.1 


Total error: 


9.0 
-8.8 


11.0 
-10.9 


12.0 
-12.1 


10.2 
-10.3 


11.1 
-10.5 


12.9 
-12.4 


±12.0 



Table 3: Systematic uncertainties (in percent) on the partial branching fractions for 
the various phase space regions, for the Babar and Belle recoil analyses. 



3 \Vub\ determinations 

The value of \Vub\ is related to the measured partial branching fractions by 



IK 



ub 



AB{B ^ XJu) 

\1 TB-AT 



theory 



a: 



where the B — > X„^z/ width according to the applied cuts, AT theory, and its un- 
certainty are determined by four theoretical calculations [SI 13 1101 IH] • Theoretical 
uncertainties can be divided in parametric terms, due to uncertainties on heavy quark 
parameters and as, and non-parametric contributions due, for instance, to higher or- 
der terms in the heavy quark expansion, weak annihilation, leading and subleading 
shape functions, renormalization scale. 

The procedure for performing the averages is documented in [7]. The average B 
lifetime used is 1.578 ps. The input values for the heavy quark parameters have been 
determined by a global fit in the kinetic scheme, translated to the scheme needed by 



Phase Space Region 


- 1) 


90% C.L. on -fwA/T 


Rl 


-0.020±0.066±0.003 


[-0.13,0.09] 


R2 


0.071±0.117±0.011 


[-0.12,0.26] 


R4 


0.042±0.066±0.009 


[-0.07,0.15] 


R5 


0.109±0.157±0.019 


[-0.15,0.37] 



Table 4: Results for {R 
under study by Babar. 



/° — 1) and limits on Jwa/'^ for the various kinematic regions 
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Kin. Expt. BLNP DGE GGOU ADFR 

region g 9^ ^ [11^ 

Endpoint analyses 

El Babar 4.35 ± 0.25l° i^J 4.15 ±0.28^011 4.17 ± 0.24t°-20 3.98 ±0.27toit 

E2 Belle 4.81 ± 0.45l^;g 4.66 ± 0.43l|5;25 4.65 ± 0.43lo;3^ 4.53 ± 0.42l° J^ 

E3 CLEO 4.00 ± 0.47l[5:^4 3.70 ± 0.43l[!:^^ 3.81±0.Ut°oil 3.47 ± 0.4ll[J J^ 

E4 Babar 4.48 ± 0.30l|i:p 4.15 ± 0.28±|j|ij n.a. 3.87 ± 0.26±g|4 

Sim.ann. Belle 4.39 ± 0.461^;^^ 4.30 ±0.45+'^- i.^ 4.24 ± 0.451^;^;^ 3.94 ± 0.4ll^ j^ 
Recoil Analyses 

Rl Babar 4.03 ± 0.19l^-26 4.23 ± 0.20t|^ ig 3.96 ± 0.18t|]-27 3.86 ± O.lSl^ ^s 

R2 Babar 3.92 ± 0.22l°;25 4.04 ± 0.22l°;26 3.84 ± 0.2ltg;^^ 3.78 ± 0.2ll[JJi^ 

R3 Babar 3.90 ± 0.24l°;28 3.93 ± 0.24l°;^^ 3.64 ± 0.22^° j° 3.60 ± 0.22l°;23 

R4 Babar 4.22 ± 0.22l°-^o 4.10 ± 0.22j:°-23 4.07 ± 0.22+°-^4 3.78 ± 0.20j;°'^i^ 

R5 Babar 4.27 ± 0.24l°-^3 4.34 ± 0.24l°-^^ 4.29 ± 0.24^°-}^ 4.34 ± 0.24;^'^^ 

R6 Babar 4.22 ± 0.27l|^;2? 4.27 ± 0.27l|^;i^ 4.21 ± 0.27t;^;ig 4.28 ± 0.27l^;25 

R5 Belle 4.45 ± 0.27±[i|i 4.53 ± 0.27t[j:|^ 4.47 ± 0.27±g;i5 4.55 ± 0.30±|j:27 

Average 4.30 ± 0.161^-23 4.37 ± 0.15lol6 4.30 ± O.ielJJ^o 4.05 ± 0.13l^-2t 

X^/d.o.f. (CL) 12.2/11 (0.36) 7.52/11 (0.76) 12.2/10 (0.27) 28.2/11 (0.003) 

Table 5: Results for \Vub\ obtained with four theoretical calculations. The uncer- 
tainties are experimental [i.e. sum of statistical and experimental systematical) and 
theoretical, respectively. The phase space regions are defined in Tables [T] and [2J 



each model, where both b — )■ ciu and b ^ s'y moments are used, giving mh{kin) = 
4.591 ± 0.031 GeV, fiUkin) = 0.454 ± 0.038 GeV^, and a correlation of -40.5%. 

Preliminary averages from the Heavy Flavour Working Group are shown in Ta- 
ble [5], for the four available theoretical calculations. All methods give consistent 
results and comparable uncertainties. 

A recent NNLO calculation fL2\ of the leading term in the partial rates gives a 
surprising change with respect to the NLO calculation used in the BLNP method, 
therefore suggesting an underestimate of the theoretical uncertainty due to the renor- 
malization matching scale. The corresponding change in \Vub\ is of the order of 8%. 
Similar estimates for the other methods are not available yet. This effect has not 
been taken into account in the current \Vub\ determinations. 



4 Conclusion 

Several \ Vub\ determinations from measurements of inclusive spectra on the full data 
sets collected at the B Factories are now available. In the most inclusive kinematic 
regions, the dominant error is due to uncertainties in signal modeling, which directly 
propagate in the selection efficiency. This uncertainty can be reduced by applying 
stricter kinematic cuts, but theoretical uncertainties correspondingly increase. Back- 
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ground modeling dominates endpoint analyses. Statistical and systematic uncertain- 
ties on the partial rates are comparable. Contributions due to weak annihilation are 
now being addressed with data; however, the statistical sensitivity is still beyond 
theoretical expectations. 

The uncertainty on inclusive \Vub\ determinations is at the 6% level, dominated 
by parametric errors (4% from a 40 MeV uncertainty on the b quark mass). Determi- 
nations by using the four available calculations are consistent between each other; all 
methods give comparable theory uncertainties, and the spread among calculations is 
comparable to theory errors. Recent results on NNLO calculations might hint to an 
underestimate of non-parametric theory uncertainties, which needs to be clarified. 
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